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D ' It is shown that monotonic depletion of the magnetic field pressure in a 
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quaking neutron star undergoing Lorentz-force-driven torsional seismic vibra- 
tions about axis of its dipole magnetic moment is accompanied by the loss of 
vibration energy of the star that causes its vibration period to lengthen at a rate 



hH ■ proportional to the rate of magnetic field decay. Highlighted is the magnetic- 



field-decay induced conversion of the energy of differentially rotational Alfven 



Qh' vibrations into the energy of oscillating magneto-dipole radiation. A set of rep- 

O ■ resentative examples of magnetic field decay illustrating the vibration energy 



powered emission with elongating periods produced by quaking neutron star are 
considered and discussed in the context of theory of magnetars. 
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Most, if not all, reported up to now computations of frequency spectra of poloidal and 
O ' toroidal Alfven vibration modes in pulsars and magnetars rest on tacitly adopted assumption 

^-^ ■ about constant-in-time magnetic field in which a perfectly conducting neutron star matter 

undergoes Lorentz-force-driven oscillations (e.g. Carroll et al. 1986, Bastrukov et al. 1997, 
Lj - 1999, Lee 2007, 2008, Shaisultanov, Eichler 2009, Sotani et al. 2009; see, also, references 

Vh ' therein). A special place in the study of the above Alfven modes of pure shear magneto- 

mechanical seismic vibrations (a-modes) occupies a homogeneous model of a solid star with 
the uniform density p and frozen-in poloidal static magnetic field of both homogeneous and 
inhomogeneous internal and dipolar external configuration. This field can be conveniently 
represented in the form 

B(r) = B b(r), b(r) = [6,.(r) ^ 0, ^^(r) ^ 0, 6^(r) = 0], B = constant, (1) 
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where B is the field intensity [in Gauss] and b(r) stands for tlie dimensionless vector-function 
of spatial distribution of the field. In particular, the above form of B(r) has been utilized 
in recent works (Bastrukov et al. 2009a, 2009b, 2009c) in which the discrete spectra of 
frequencies of node-free torsional Alfven oscillations has been computed in analytic form on 
the basis of equations of linear magneto-solid mechanics (Bastrukov et al. 2009a) 

pu = -[6}xB], V-u = 0, (2) 

c 

5j = -^[VxdB], 5B = Vx[uxB], V-B = 0. (3) 

4:77 

Here u = u(r, t) is the field of material displacement, the fundamental dynamical vari- 
able of solid-mechanical theory of elastically deformable (non-flowing) material continua. 
These equations describe Lorentz-force-driven non-compressional vibrations of a perfectly 
conducting elastic matter of a solid star with Ampere form of fluctuating current density 
5j. Equation for 5B describing coupling between fluctuating field of material displacements 
u and background magnetic field B pervading stellar material is the mathematical form of 
Alfven theorem about frozen-in lines of magnetic field in perfectly conducting matteiLj. As 
for the general asteroseismology of compact objects is concerned, the above equations seems 
to be appropriate not only for neutron stars but also white dwarfs (Molodtsova et al. 2010) 
and quark stars (Heyvaerts et al. 2009). The superdense material of these latter yet hypo- 
thetical compact stars is too expected to be in solid state (Xu 2003, 2009), that is, in the 
aggregate state possessing mechanical property of shear elasticity. 



^Equation ([2]) can be represented in the following equivalent tensor form (Bastrukov et al. 2009a), 
Ui = ^kTik, where stress tensor Tik — {l/4:w)[Bi5Bk+Bk5Bi — BjSBjSik] is the Maxwellian tensor of magnetic 
field stresses. This last equation for m is identical in appearance to canonical equation of solid-mechanics 
Ui — "^kOik where CTifc = 2/iMifc-|-[K— (2/3)/i] u^j^ifc is the Hookean stress tensor and u^fc = (l/2)[Vi Ufc+Vfc Mi] 
is the tensor of shear deformations in an isotropic elastic continuous matter with shear modulus /i and 
bulk modulus n. The adopted for Alfven vibrations of compact solid stars terminology [which is not new 
of course, e.g.. Moon (1984)], namely, magneto-solid mechanics and solid-magnetics, is used as a solid- 
mechanical counterpart of well-known terms like magneto-fiuid mechanics, magnetohydrodynamics (MHD) 
and hydromagnetics, pSir = {l/c)[5j x B], <5j = (c/47r)[V x (5B], SB = V x [Sv x B], describing fiowing 
magento- active plasma in terms of the velocity of fiuctuating flow 5v = u and fluctuating magnetic field S'B. 
The MHD approach is normally utilized in astrophysics of main-sequence (MS) liquid stars. The advantage 
and virtue of the Rayleigh energy method of computing frequencies of node- free Alfven vibrations (Bastrukov 
et al. 2009a, 2009b, 2009c, Molodtsova et al. 2010) is that it can too be applied to the above MHD equations 
underlying the study of a-modes in liquid MS-stars such, for instance, as rapidly oscillating Ap (roAp) stars, 
chemically peculiar magnetic stars exhibiting high-frequency oscillations that have been and still remain the 
subject of extensive investigation (e.g., Ledoux and Walraven 1958, Rincon and Rieutord 2003). 
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Inserting ([2]) in ([3]) the former equation of solid-magnetics takes the form 

pu(r,t) = i-[V X [V X [u(r,t) x B(r)]]] x B(r). (4) 

The obtained equation serves as a basis of our further analysis. The studied in above works 
vibration regime is of some interest in that the rate of differentially rotational material 
displacements 

ti(r, t) = [a;(r, t) x r], a;(r, t) = [Vx(r)] a{t), (5) 

V\{v) = 0, x{r) = fiPiicos9), fe{r) = A,/ (6) 

has one and the same form as in torsion elastic mode of node-free vibrations under the action 
of Hooke's force of mechanical shear stressecl (Bastrukov et al. 1999, 2007, 2008). Hereafter 
P^(cos0) stands for Legendre polynomial of degree i specifying the overtone of a- mode. 
The time-dependent amplitude a(t) describes temporal evolution of above vibrations; the 
governing equation for a{t) is obtained form equation (jl]). The prime purpose of above works 
was to get some insight into difference between spectra of discrete frequencies of toroidal a- 
modes in neutron star models having one and the same mass M and radius R, but different 
shapes of constant-in-time poloidal magnetic fields. By use of the energy method, it was 
found that each specific form of spatial configuration of static magnetic field about axis of 
which the neutron star matter undergoes nodeless torsional oscillations is uniquely reflected 
in the discrete frequency spectra by form of dependence of frequency upon overtone i of 
nodeless vibration. 

The assumption about constant in time undisturbed magnetic field means that the 
internal magnetic field pressure, Pb, the velocity va of Alfven wave in the star bulk (e.g. 



^It may be worth noting that in standing- wave regime of elastic- force-driven differentially rotational 
vibrations, the toroidal field of material displacements is described one and the same equation ([5]), but 
function x = Aiji{kr)Pi {cos 9), is the solution of Helmholtz equation V^x -I- k^x — where fc^ — uj^/c^ and 
c^ = /i/p is the speed of shear wave of material displacements in solid stellar matter of density p and whose 
shear modulus is /i. The discrete frequency spectrum uj{i) = 2-kv[I!.) of elastic torsional vibration of multipole 
degree (. is computed from from boundary condition 71^(7^. = resulting in dispersion equation, dji{z)/dz = 
ji{z)l z, where z = kR. The characteristic feature of the long-wave regime of node-free vibrations, z << 1, 
with the toroidal field of material displacements ©-([S]) is the absence of dipole overtone in the spectrum of 
discrete frequencies (Bastrukov 2007a, 2007b): J^e(-^) = Ve[{2^ + 3)(^ — 1)]"^^^, with basic frequency of shear 
elastic vibrations, Ve — ct/R, where R stands for the star radius. In the non-studied before regime of short- 
wave length limit of torsional elastic vibrations (the limit z >> 1 taken in the exact dispersion equation 
for dipole overtone dji{z)/dz = ji[z)/z) leads to the spectrum of standing- wave dipole vibrations of the 
following analytic form Vf,{l = 1) = Vf.{p + 1) , where p = 0,l,2,3... is the node number of short- wavelength 
vibration regime. 



-4- 



10^ 



10' 



£ 10'' 

< -1 
> 10 



10" 



10" 



— 1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I ■: 



M=1.2M ,R=10km 

sun' 

M=0.8M , R=20 km 

sun' 



J I I I I I I I I I I I 



M=1.2M3^„,R=10km 

M=0.8M , R=20 km 

sun' 




— 1 I I I I I I I I I I I 



10 



12 



10 



13 



10 



14 



10 



15 



10 



16 



B 



Fig. 1. — The basic frequency and period of global Alfven oscillations, equations (IE]), as 
functions of magnetic field intensity in the neutron star models with indicated mass and 
radius. 
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and, hence, the frequency va = wa/Stt (where oja = va/R) and the period Pa 
global Alfven oscillations 
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(8) 



remain constant in the process of vibrations whose amplitude a{t) subjects to standard 
equation of undamped harmonic oscillator (e.g., Bastrukov et al 2009a). The allow for 
viscosity of stellar material leads to exponential damping of amplitude, but the frequency 
ua and, hence, the period Pa preserve one and the same values as in the case of non- viscous 
vibrations (Bastrukov et al. 2009c). In Fig.l these latter quantities are plotted as functions of 
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intensity B of undisturbed poloidal magnetic field in the neutron star models with indicated 
mass M and radius R. The practical usefulness of chosen logarithmic scale in this figure is 
that it shows absolute vales of z/^ and P4 for global vibrations of typical in mass and radius 
neutron stars. 

In this work, continuing a line of investigations reported in (Bastrukov et al. 2009a, 
2009b, 2009c, Molodtsova et al. 2010), we relax the assumption about constant-in-time mag- 
netic field and examine the impact of its decay on the vibration energy and period. In section 
2, a mathematical background of quaking neutron star model is briefly outlined with empha- 
sis on the loss of vibration energy caused by depletion of internal magnetic field pressure and 
resulting vibration-energy powered magneto-dipole radiation. The decreasing of magnetic 
field pressure in the star is presumed to be caused by coupling between vibrating star and 
outgoing material which is expelled by quake, but mechanisms of star-envelope interaction 
resulting in the decay of magnetic field, during the time of vibrational relaxation, are not 
considered. It is shown that physically meaningful inferences regarding radiative activity of 
quaking neutron star can be made even when detailed mechanisms of depletion of magnetic 
field pressure in the process of vibrations triggered by quakes are not exactly known. In sec- 
tion 3, a set of representative examples illustrating the newly disclosed effects is considered. 
The summary of developed theory with emphasis on its relevance to electromagnetic activity 
of quaking magnetars is presented in section 4. 



2. Vibration-energy powered magneto-dipole emission of neutron star 

In what follows we focus on dependence on time of the intensity of poloidal magnetic 
field, B = B{t), about axis of which the star undergoes global torsional vibrations and 
confine our consideration to the model of uniform internal field 

B(r,t) = fi(t)b(r), [br = cosO, bg = -sinO, b^ = 0]. (9) 

On account of this the equation of solid- magnetics, (jl]), takes the form 

pii(r,t) = ^[V X [V X [u(r,t) x b(r,t)]]] x b(r). (10) 

Inserting here the following separable form of fluctuating material displacements 

u(r,t) = a(r)a(t) (11) 

we obtain 

{pa(r)}«(t) = 2P^(t){[V X [V X [a(r) x b(r)]]] x b(r)}a(t), Psit) = ^.(12) 
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Scalar product of flT2|) with the time-independent field of instantaneous displacements a(r) 
followed by integration over the star volume leads to equation for amplitude a{t) having the 
form of equation of oscillator with depending on time spring constant 

Ma{t) + }C{t)a{t)=0, M=pme, IC{t) = 2PBit) ke, (13) 

mi= a(r) ■ a(r) dV, a = AV x [r/P^(cos^)], (14) 

ke = /"a(r) • [b(r) x [V x [V x [a(r) x b(r)]]]] dV. (15) 

The solution of equation of non- isochronal (non-uniform in duration) and non-stationary 
vibrations with time-dependent frequency [a(t) -|- u^(t)a(t) = where w^(t) = /C(t)/A^] is 
non-trivial and fairly formidable task (e.g., Vakman and Vainshtein 1977). But solution of 
such an equation, however, is not a prime purpose of this work. The main subject is the 
impact of depletion of magnetic-field-pressure on the total energy of Alfven vibrations 

E4t) = ^ + ^^ (16) 

and the discrete spectrum of frequency of the toroidal a-mode 

ujm=u:l{t)sl c.i(t) = ^, s^, = ^R\ (17) 



u]{t) = B\t)Kl k] = -^, s] 



AixpR^ ' 



{e - 1 
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It is to be stated clearly from the onset that it is not our goal here to speculate about possible 
mechanisms of neutron star demagnetization and advocate conceivable laws of magnetic field 
decay. The main purpose is to get some insight into the effect of arbitrary law of magnetic 
field decay in quaking neutron star on period of Lorentz-force-driven seismic vibrations and 
radiative activity of the star brought about by such vibrations. In the reminder of the paper 
we remove index i in the expression for discrete spectrum of frequency of toroidal Alfven 
mode putting uj{t) = ue{t). It seems worth noting that first computation of discrete spectra 
of frequencies of toroidal Alfven stellar vibrations in standing wave- regime, has been reported 
by by Chandrasekhar (1956). The extensive review of other earlier computations of discrete 
frequency spectra of a-modes, ui = ojaSi, is given in well-known review of Ledoux and 
Walraven (1958). In present work, we follow the line of the energy method developed in works 
(Bastrukov et al 2009a, 2009b, Molodtsova et al 2010) for computing discrete frequencies of 
node-free Alfven vibrations and preserve one and the same notations for frequency spectra 
as in this latter canonical paper on variable stars. 



by 
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2.1. Magnetic-field-decay induced loss of vibration energy 

The total energy stored in quake-induced Alfven seismic vibrations of the star is given 

E,(() = ^1^ + ^(M)f!!W , ^B(ty,=.\BmiM. (19) 



Perhaps most striking consequence of the magnetic-field-pressure depletion during the post- 
quake vibrational relaxation of neutron star is that it leads to the loss of vibration energy 
at a rate proportional to the rate of magnetic field decay 



= Ate!W M^^ ^ ^,VWBW^. (21) 

2 dB at at 

It is worth emphasizing at this point that the magnetic-field-decay induced loss of vibration 
energy is substantially different from the vibration energy dissipation caused by shear vis- 
cosity of matter resulting in heating of stellar material (Mestel 1999). As was noted, the 
characteristic feature of this latter mechanism of vibration energy conversion into the heat 
(i.e., into the energy of non-coherent electromagnetic emission responsible for the formation 
of photosphere of the star) is that the frequency and, hence, period of vibrations are the same 
as in the case of viscous- free vibrations (Bastrukov et al. 2009c). However, it is no longer so 
in the case under consideration. It follows from above that depletion of magnetic field pres- 
sure resulting in the loss of total energy of Alfven vibrations of the star causes its vibration 
period to lengthen at a rate proportional to the rate of magnetic field decay. The process of 
magnetic field decay depends on many factors of demagnetization of neutron star matter and 
likely primarily on the star-environment communication. Also, it seems quite likely that, 
contrary to viscous dissipation, the loss of vibration energy due to decay of magnetic field 
must be accompanied by coherent (non-thermal) electromagnetic radiation. Adhering to this 
supposition in the next sections we consider several representative scenarios of magnetic field 
decay in quaking neutron stars. In so doing, special emphasis is made on conversion of the 
energy of Lorentz-force-driven seismic vibrations into the energy of magneto-dipole emission 
whose flux oscillates with frequency of torsional Alfven magneto-mechanical vibrations of 
the star. 



2.2. Conversion of vibration energy into power of magneto-dipole radiation 

The point of departure in the study of vibration-energy powered magneto-dipole emis- 
sion of the star (whose radiation power, P, is given by Larmor's formula) is the equation 

'-Ml.r. P.A.K,,. (22) 

Consider a model of quaking neutron star whose torsional magneto-mechanical oscillations 
are accompanied by fluctuations of total magnetic moment preserving its initial (in seismi- 
cally quiescent state) direction: fi = fin = constant. The total magnetic dipole moment 
should execute oscillations with frequency uj{t) equal to that for magneto-mechanical vibra- 
tions of stellar matter which are described by equation for a{t). This means that 6fi{t) and 
a{t) must be subjected to equations of similar form, namely 

6il{t)+uj\t)Sn{t)=0, (23) 

a{t) + uj^{t)a{t) = 0, tu^it) = B'^{t)K^. (24) 

It is easy to see that equations (|27|) and fl28|) can be reconciled if 

6fj,{t) = ina{t). (25) 

Then, from ( 123|) . it follows 6fl = —iuo'^^a. Given this and equating 

^^ = A^«:V(t)i?(t)^, (26) 

dt ^ ' ^ ' dt ^ ' 

with 

V = -^^fi'K'B\t)a\t) (27) 

we arrive at the equation of time evolution of magnetic field 

dBit) ., , 211^^^ 

which yields the following law of its decay 

5(t) = ^£E=, ^-1= 27^^(0). (29) 

The lifetime of magnetic field r is regarded as a parameter whose value is established from 
above given relations between the period P and its time derivative P which are taken from 
observations. And knowing from observations P{t) and P{t) and estimating r one can get 
information about the strength of magnetic field in the star and the magnitude of its total 
magnetic moment in undisturbed statqj. 



•^Understandably that in the model of vibration-energy powered magneto-dipole emission under consid- 
eration, the equation of magnetic field evolution is obtained in similar fashion as equation for the angular 



-9- 




>. 0.2 

X 



-• 
-0.2. 



- P(x)/P(0) 

- AP(x)/P(0) 



P(x)/P(0) 

AP(x)/P(0) 







0.2 



0.4 0.6 

x = t/x 



0.8 



Fig. 2. — The upper panel illustrates the magnetic-field-decay induced elongation of the 
period P{x) and difference between periods AP{x) of toroidal a- mode, defined in the text, 
plotted as functions of x = t/r. The lower panel shows that in this model of magnetic 
field decay the time derivatives of period P{x) and difference between periods AP{x) are 
decreasing functions of x. 



2.3. Lengthening of vibration period 

The immediate consequence of above line of argument is the magnetic-field-decay in- 
duced lengthening of vibration period 



Pit) 



2tt 



C 



c = — 

uj{t) B{t)' k' 



<0. 



The rate of period elongation is given by 

dP{t) _ C dB{t) dB{t) 

dt ~ ^ B^{t) dt ' dt 

Combining these equations we obtain the following general relations 

P{t)B{t) = constant, 

Pjt) _ Bit) 
Pit) Bit) 



(30) 

(31) 

(32) 
(33) 



velocity 51 does in the standard model of rotation-energy powered emission of neutron star. However, in the 
model of quaking neutron star vibrating in toroidal a-mode, the elongation of period is attributed to mag- 
netic field decay, whereas in canonical Pacini-Gold model of radio-pulsar the period lengthening is attributed 
to the spin-down of rotating neutron star. 
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which are generic to solely Alfven vibration mode and independent of specific form of the 
magnetic field decay law. The difference between periods evaluated at successive moments 
of time ti = and t2 = t is given by 



AP(t) = P(t) - P(0) = -P(0) 



1 



m 

B{t) 



> 0, P(t) < P(0). 



(34) 



The practical usefulness of these general relations is that they can be used as a guide in search 
for fingerprints of Alfven seismic vibrations in data on oscillating emission from quaking 
neutron star. 



In the model under consideration 

Pit) 
P(0) 
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wr 
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wr 



Bit) 



P(0) 



P(0) = B{t = 0). 



and the time evolution of vibrational period and its derivative is described by 



P(t) = P(0) [1 + (t/r)]i/2, AP(t) = P(0) 1 - ^/iTtJ^ 



Pit) 



P(0) 



2r[l + (t/r)]V2' 



AP(t) = P{t) - P(0) 



It follows that lifetime r is determined by 



^ P(0) 



P(0) 



2r 
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2r 
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VT+t/'^i 



(35) 
(36) 

(37) 
(38) 

(39) 



P2(0) 
P(t) P(t) = _-^^ = constant 



and the ratio P to P is given by 



Pit) 2r 



2r 



[1 + itlr)]-\ 



(40) 



(41) 



For the sake of illustration of general trends in P and P, in Fig. 2 these quantities are plotted 
as functions of the fraction time x = t/r which is ranged in the interval < a; < 1. In Fig. 
3, we plot P(t)/P(0) versus P(t)/P(0). 
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Fig. 3. — The interrelation between time derivative of period P{x) versus period P{x), both 
normahzed to period at initial moment of time -P(O), computed as functions of x = t/r. 

2.4. Time evolution of the vibration amplitude 

It is interesting to note that considered model permits exact analytic solution of equation 
for vibration amplitude a{t) which is convenient to represented as 



a{t)+uj'^{t)a{t) =0, uj^{t) 



u\0) 



l + t/r 
The procedure is as follows. Let us introduce new variable 



, u{0) = uai^ = constant. 



s = l + t/r, -^ t=(s-l)r 
dt 



T 



a{t) 



dt 



a'is) 



da{s) 



ds 



so that 



a(t) = ^ = ^^ = ^, a{t) 



a" is) 



T 



dt ds dt 

In terms of a{s), equation (H2]) takes the form 

sa"{s) + I3^a{s) = 0, jS^ = u^{0)t'^ = constant. 

This equation permits exact analytic solution (e.g. Polyanin and Zaitsev, 2004) 

a(s) = s'/^{Ci Ji{2(3s^'^) + C2Y^{2(3s^'^)} 

where Ji(2/3s^/^) and Yi{2(5s^/'^) are Bessel functions (Abramowitz and Stegun 1972) 

J^{z) = - [ cos(zsin^ - 9)d9, Yi{z) = - [ sin(2sin^ - 9)d9, z = 2l3s^/^ 

T^ Jo T^ Jo 



(42) 

(43) 
(44) 

(45) 
(46) 
(47) 

(48) 
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The arbitrary constants Ci and C2 can be eliminated from two conditions 

a{t = 0) = tto, a{t = r) = 0. (49) 

where zero-point amplitude 

(50) 



a, 



2 _ 2Ea{0) _ 2Ea{0) ^2(q^ _ K{0) 



is related to the average energy Ea{0) stored in torsional Alfven vibrations at initial (before 
magnetic field decay) moment of time t = 0. So, knowing from observations the amplitude 
of vibrations Oq, the average energy £'^(0) is established from last equationc|. 
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Fig. 4. — Vibration amplitude a{t) computed at fixed C and (3 and different values of rj. 
This shows that variation of this parameter is manifested in change of magnitude of | a | , but 
period elongation is not changed. 



As a result, the general solution of ( l42l) can be represented in the form 

a{t) = [1 + (t/r]i/2C{ Ji(2/3 [1 + (t/r)]i/2) - r^n(2/3 [1 + {t/ry/'])}, 
JMr)) 



r] 



YMT)y 



C = ao[Mz{0))-r]Y,{z{0)]-\ 



(51) 
(52) 



^Before magnetic field decay, the star oscillates in the harmonic in time regime, that is, with amplitude 
a = Q;ocos(a;(0)t), so that a^ = (l/2)aQ and d^ = {l/2)uj^{0)aQ. The average energy stored in such 
oscillations is given by EAiO) = (l/2)Md2 + {l/2)K{0)a^ = (l/2)Mw2(0)ag = {l/2)K{0)al from which 
equation P^ is obtained. 
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The vibration period lengthening in the process of vibrations is illustrated in Fig.4 and 
Fig.5, where we plot a{t), equation fjSTj) . at different values of parameters (3 and r] pointed 
out in the figures. Fig.4 clearly shows that rj is the parameter regulating magnitude of 
vibration amplitude, the larger rj, the higher is amplitude. However this parameter does not 
affect the rate of period lengthening. As it is clearly seen from Fig.5 both the elongation rate 
of vibration period and magnitude of vibration amplitude are highly sensitive to parameter 
/3. 
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Fig. 5. — Vibration amplitude a(t) computed at fixed C and rj and different values of /3. 



3. Representative examples 



The considered law of magnetic field decay cannot be, of course, regarded as universal 
because it reflects a quite concrete line of argument regarding the fluctuations of magnetic 
moment of the star. The interrelation between quake-induced oscillations of total magnetic 
moment of the star 6fi{t) and the amplitude a{t) of its seismic magneto- mechanical oscilla- 
tions can be consistently interpreted with the aid of the function of dipole demagnetization 
f{B{t)) which is defined by the following condition of self-consistency in a of right and left 
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hand sides of equation fl22|) . namely 

6fl{t) = ii{B{t))a{t) (53) 

The vector-function of dipole demagnetization f[B{t)) depending on decaying magnetic field 
reflects temporal changes of electromagnetic properties of neutron star matter as well as 
evolution of magnetic-field-promoted coupling between neutron star and its environment. 
This means that specific form of this phenomenological function should be motivated by 
heuristic arguments taking into account these factors. With this form of Sil{t), equation 
(!22|) is transformed to magnetic field decay of the form 

dB{t) f'(5(t)) 2 

In the above considered case this function is given by This implies that the function of dipole 
demagnetization f{B{t)) has the form 

f (5(t)) = (3 B{t) B(t), /3 = kV = constant. (55) 

As is shown below, the practical usefulness of function of dipole demagnetization, f (i?(t)), 
consists in that it provides economic way of studying a vast variety of heuristically motivated 
laws of magnetic field decay, B = B{t), whose inferences can ultimately be tested by obser- 
vations. In the next section, with no discussion of any specific physical mechanism which 
could be responsible for magnetic field decay, we consider a set of representative examples of 
demagnetization function i{B{t)) some of which lead to the magnetic field decay laws that 
have been regarded before, though in a somewhat different context (e.g. Sang, Chanmugam 
1987, Srinivasan et al. 1990, Goldreich, Reisenegger 1992, Urpin et al. 1994, Wang 1997, 
Livio et al. 1998, Geppert et. al 2001, Bisnovatyi-Kogan, 2002). In this connection it worth 
stressing that the model under consideration deals with magnetic field decay in the course 
of vibrations triggered by starquake, not with long-term secular decay which has been the 
subject of these latter investigations (see also references therein). In the case in question, 
a fairly rapid decay of magnetic field during the time of post-quake vibrational relaxation 
of the star is though of as caused, to a large extent, by coupling of the vibrating star with 
material expelled by quake. In other words, escaping material removing a part of magnetic 
flux density from the star is considered to be a most plausible reason of depletion of inter- 
nal magnetic field pressure in the star. While physical processes responsible for magnetic 
field decay remain uncertain, practical significance of above, admittedly heuristic, line of 
argument is that it leads to meaningful conclusion (regarding elongation of periods of os- 
cillating magneto-dipole emission) even when detailed mechanisms of depletion of magnetic 
field pressure in the course of vibrations are not exactly known. 
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Fig. 6. — The upper panel - the period P[x) and difference between periods AP(x) of 
toroidal a-mode as functions of x = t/r. The middle panel shows that in this model of 
magnetic field decay the time derivatives of period P[x) and difference between periods 
AP(a;) are decreasing functions of x. The lower panel, the caused by field decay rate of 
the a-modes vibration period lengthening, normalized to period of initial moment of time, 
caused by magnetic field decay with time evolution of considered in the text representative 
examples i = \ and 2. 



1. As a first representative example, consider a model of quaking neutron star in which 
function of dipole demagnetization is given by 



f(5(t)) = /3B(t). 

This leads to exponential decay of magnetic filed 
dB{t) B{t) 



B{t) = B{0) e-'l\ T-^ = T]P^ 



(56) 



(57) 



dt T 

where -B(O) is the intensity of magnetic field before quake. The period and its derivative are 
given by 



P(t) = P(0) e*/^ P{t) 



P(0) 



e*/^ 



T 



(5^ 
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Fig. 7. — The rate of period versus period of toroidal a-mode computed in 1-st and 2-nd 
representative models. 

The difference between periods, evaluated at the moment of time t2 = t and ti = is given 
by 



AP(t) = P{t) - P(0) = P(0) e*/"^(l - e"*/") 



(59) 



In the initial stage of decay, when (t/r) << 1, the last equation is reduced to AP(t) 
P(0) [t/r) [1 + [t/r)]. The parameter of magnetic field lifetime is given by 



r 



p{ty 



(60) 



In Fig. 6 the above period and its derivative are plotted as functions of fraction time x = t/r. 
2. Consider a neutron star model whose function of demagnetization is given by 

(61) 



f(P(t))=/3v/p(tyB(t). 

The resultant equation of magnetic field evolution reads 

dB{t) 



dt 



-^B\t), 7 = ^7/3' 



(62) 



and thus implying that rate of magnetic field decay is proportional to the density of magnetic 
field energy stored in the star. In this case one has 



B{t) = P(0) 1 + 



t 



The elongation of vibration period and a rate of its lengthening are given by 
P(t) = P(0) (l + - J , P(t) = —^ = constant. 



(63) 



(64) 
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The difference between periods is given by AP(t) = P(0)(t/r). The decay time is eliminated 
from the ratio |j 



Pit) 



t 



Pit) ^'v^rj ^ p{oy 



1 Pio) 



(65) 



The inferences of this model regarding period and its derivative as functions of time are 
presented in Fig. 7. 




Fig. 8. — The rate of frequency z>(a;) normalized to i/(0) as a function of x = t/r for 
logarithmic law of magnetic filed decay. 



3. Finally, consider a model with quite sophisticated function of dipole demagnetization 

(66) 



i{B{t)) 



BiO) B(t) 



r + t,/;^Bit) 



which is interesting in that it leads to equation 

dBjt) _ BjO) / ^ t 

dt T \ T 



BiO) 

T + t 



^Similar analysis can be performed for the demagnetization function of the form 



which leads to 



{(Bit)) = /3v/S'"-i(t)B(i), m = 3, 4, 5... 
dB{t) 



dt 
B{t) = 



= -7i3"(t), 7 = ^/3', 



= 7(m-l)B™-i(0). 



(67) 
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yielding a fairly non-trivial logarithmic law of magnetic field decay 



Bit) = 5(0) 



1 - In 1 + 



(68) 



The corresponding frequency of toroidal a-mode as a function of time is plotted in Fig. 8 
and it may be interesting to note that similar shape of z>(t)/z/(0) exhibit data on post-glitch 
emission of PSR J1846-0258 (Livingstone et al. 2010). 



In this model 



P{t)- 


= P{0) 




-1 
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m- 


dPit) 
dt 


_P(0) 

T + t 


1-ln 


f t\\ 

1 + - 

V ^/. 



no)" 



(69) 
(70) 
(71) 



These examples show that while the period elongation is the common effect of depletion of 
magnetic field pressure, the interrelations between periods and its derivatives substantially 
depend on specific form of the magnetic field decay law. 



4. Summary 

The main purpose of this work was to examine inferences of asteroseismic model of 
neutron star undergoing quake-induced torsional Alfven vibrations about axis of its dipole 
magnetic moment which are accompanied by monotonic depletion of internal magnetic field 
pressure. It is shown this latter process results in the loss of vibration energy and changes 
the character of temporal evolution of vibrations from harmonic in time vibrations (with 
constant frequency) in static, time-independent, magnetic field to non-isochronal vibrations 
with monotonically decreasing frequency in decaying magnetic field. The presented analysis 
suggests that the loss of vibration energy caused by depletion of magnetic field pressure 
should be accompanied by coherent (non-thermal) oscillating magneto-dipole radiation (with 
frequency of luminosity oscillations equal to the frequency of Alfven seismic vibrations of 
the star). In this context it is appropriate to remind a seminal work of Hoyle, Narlikar and 
Wheeler (1964) in which it has been pointed out for the first time that vibrating neutron 
star should operate like Hertzian magnetic dipole deriving radiative power of magneto-dipole 
emission from the energy of magneto- mechanical vibrations (see, also, Pacini 2008). What is 
newly disclosed here is that conversion of vibration energy into the energy of magneto-dipole 
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electromagnetic radiation can be realized when, and only when, torsional Alfven vibrations 
are accompanied by decay of magnetic field. 

All above suggests that considered for the first time theory of vibration-energy powered 
emission of neutron star is relevant to electromagnetic activity of magnetars - neutron stars 
endowed with magnetic field of extremely high intensity the radiative activity of which is 
ultimately related to the magnetic field decay (Duncan & Thompson 1992, Thompson & 
Duncan 1995). This subclass of highly magnetized compact objects is commonly associated 
with soft gamma repeaters and anomalous X-ray pulsars (e.g., Harding 1999, Kouveliotou 
1999, Woods and Thompson 2006, Mereghetti 2008) - young isolated and seismically active 
neutron stars (Chen et al. 1996, Franco et al. 2000). The magnetar quakes are exhibited 
by short-duration thermonuclear gamma-ray fiash followed by rapidly oscillating X-ray fiare 
of several-hundred-seconds duration. During this latter stage of quake-induced radiation 
of magnetar, a long-periodic (5-12 sec) modulation of brightness was observed lasting for 
about 3-4 minutes (e.g., Harding 1999, Kouveliotou 1999). Taking into account such long- 
periodic and rapidly dying variability of magnetar pulsating radiation is the feature which 
is non-typical for young neutron stars powered by energy of rotation, it has been suggested 
(Bastrukov et al. 2002) that detected long-periodic pulsating emission is powered by the 
energy of torsional magneto-elastic vibrations of magnetar triggered by quake. In view of 
key role of ultra-strong magnetic field it is quite likely that quasi-periodic oscillations (QPOs) 
of outburst flux from SGR 1806-20 and SGR 1900-1-14 discovered in observations reported in 
(Israel et al. 2005, Watts & Strohmayer 2006, Terasawa et al. 2006) are produced by torsional 
seismic vibrations predominately sustained by Lorentz force (Bastrukov et al 2009b, 2009c) 
which are accompanied, as was argued above, by monotonic decay of background magnetic 
field. If so, the predicted elongation of QPOs period of oscillating outburst emission from 
quaking magnetars should be traced in existing and future observations. 

The authors are indebted to colleagues from Physics Department and Kavli Institute 
of Astronomy and Astrophysics, Peking University, Beijing, for valuable discussions and 
suggestions. Special thanks go to Alexey Kudryashov (Saratov University, Russia) who paid 
our attention to the existing of exact solution of obtained equation for vibration amplitude. 
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